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Abstract Big-bang nucleosynthesis (BBN) theory, together with the precise WMAP cosmic 

T ) . baryon density, makes tight predictions for the abundances of the lightest elements. Deuterium 

and 4 He measurements agree well with expectations, but 7 Li observations lie a factor 3 — 4 

i-^ , below the BBN+WMAP prediction. This 4 — 5cr mismatch constitutes the cosmic "lithium 

problem," with disparate solutions possible. (1) Astrophysical systematics in the observations 

could exist but are increasingly constrained. (2) Nuclear physics experiments provide a wealth of 

well-measured cross-section data, but 7 Be destruction could be enhanced by unknown or poorly- 

C/j ■ measured resonances, such as 7 Be + 3 He — > 10 C* — > p + 9 B. (3) Physics beyond the Standard 

Model can alter the 7 Li abundance, though D and 4 He must remain unperturbed; we discuss 

such scenarios, highlighting decaying Supersymmetric particles and time-varying fundamental 

constants. Present and planned experiments could reveal which (if any) of these is the solution 

to the problem. 
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1 Introduction 

Big bang nucleosynthesis (BBN) describes the production of the lightest nuclides- 
D, 3 He, He, and 'Li-at times ~ 1 sec to ~ 3 min after the big bang. Theoreti- 
cal predictions of the light element abundances are well-understood and rest on 
the secure microphysics of nuclear cross sections and Standard Model weak and 
electromagnetic interactions dH[2]). These predictions are in broad quantitative 
agreement with measured primordial light element abundances derived from ob- 
servations in the local and high-redshift universe. This concordance represents a 
great success of the hot big bang cosmology, and makes BBN our earliest reliable 
probe of the universe. 

BBN has dramatically changed over the past decade in response to the cos- 
mological revolution sparked by the recent flood of new observations. Notably, 
WMAP measurements of the cosmic microwave background (CMB) radiation 
have precisely determined the cosmological baryon and total matter contents 
d3U3]) , while high-redshift supernova observations reveal that the universe recently 
entered a phase of accelerated expansion dUEEE]) . These and other observations 
reveal the existence and tightly determine the abundance of both dark matter 
and dark energy on cosmic scales where they dominate the mass-energy budget 
of the universe today. Yet the nature of both dark matter and dark energy re- 
mains unknown. Thus 21st-century cosmology finds itself in a peculiar state of 
"precision ignorance;" this situation is particularly exciting for particle physics 
because both dark matter and dark energy demand physics beyond the Standard 
Model. 

BBN has played a central role in the development of this new cosmology. CMB 
data now measure of the cosmic baryon density, independently of BBN, and with 
high precision. This casts BBN in a new light: a comparison of these two measures 
of cosmic baryons provides a strong new test of the basic hot big bang framework 
([9]). Moreover, this test can now be performed in a new way, using the CMB 
baryon density as an input to the BBN calculation, which outputs predictions for 
each light element abundance; these can then each be directly compared to light 
element observations (jlOh . The result is that deuterium shows spectacular agree- 
ment between BBN+CMB predictions and high-redshift observations, and 4 He 
shows good agreement. However, using the first-year WMAP data, 7 Li showed a 
discrepancy of a factor 2 — 3, representing a 2 — 3<r disagreement between observa- 
tions and theory (jlip i2 1 [T3 |ll4p . This disagreement has worsened over time, now 
standing at a factor 3 — 4 in abundance or 4 — 5a: this is the "lithium problem." 

In this paper we present an overview of the lithium problem, accessible to 
nuclear and particle physicists and astrophysicists. Broader reviews of primordial 
nucleosynthesis and its relation to cosmology and particle physics are available 

In Section [21 we trace the origin of the lithium problem, with a focus on 
the physics of BBN 7 Li production, the nature and precision of light element 
abundance measurements, and the state of light element concordance in view of 
the CMB-measured cosmic baryon density. We review possible solutions to the 
lithium problem in Section [3) (i) astrophysical systematic uncertainties in lithium 
abundances and/or their interpretation; or (ii) new or revised nuclear physics in- 
puts to the BBN calculation, in the form of increased mass-7 destruction via novel 
reaction pathways or by resonant enhancement of otherwise minor channels; or 
(iii) new physics - either particle processes beyond the Standard Model occur- 
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ring during or soon after BBN, or large changes to the cosmological framework 
used to interpret light element (and other) data. We close by summarizing the 
near-future outlook in Section [U 

2 Standard BBN in Light of WMAP: the Lithium Problem Emerges 

2.1 Standard BBN Theory 

The cosmic production of light nuclides is the result of weak and nuclear reactions 
in the context of an expanding, cooling universe. "Standard" BBN refers to 
the scenario for light element production which marries the Standard Model of 
particle physics with the "standard" (ACDM) cosmology, with: 

1. gravity governed by General Relativity 

2. a homogeneous and isotropic universe (cosmological principle) 

3. the microphysics of the Standard Model of particle physics 

4. the particle content of the Standard Model, supplemented by dark matter 
and dark energy 

Under these assumptions, the expansion of the universe is governed by the Frei- 
dmann equation 

where a(t) is the dimensionless cosmic scale factor (related to redshift z via 
1 + z = 1/a), and H = a /a is the universal expansion rate. The total cosmic 
mass-energy density p = J2 Pi sums contributions from all cosmic species i. 

By far the largest contribution to the density comes from radiation: relativistic 
species for which m <C T (with T the temperature), namely blackbody photons 
and N v = 3 species of neutrinos and antineutrinos, and e^ pairs at T > m e . 
Cosmic matter consists of nonrelativistic species with m S> T: nucleons n and p, 
and e~ at T < m e . Since /9 ra d S> p m attcr, eq. CEJ shows that radiation dominates 
cosmic dynamics during BBN. 

BBN occurs entirely in this radiation-dominated epoch, for which the energy 
density has p oc T 4 , where T oc 1/a (adiabatic cooling). This, together with 
eq. ([1]), gives t oc 1/T 2 , or 

1 MeV\ 2 



1 sec ( —^r- 1 (2) 



Light-element formation depends crucially on the relative amounts of baryons 
(nucleons) and radiation, parameterized by the baryon-to-photon ratio 

r, = ^ = 2.74 x W~ 8 n h h 2 . (3) 

Here Q^ = Pb/pcvit arid p CT - lt = 3Hq/8ttG, with Hq the present expansion rate 
(Hubble parameter). In standard BBN, rj is the only free parameter controlling 
primordial light element abundances. 

Initially, cosmic baryons are in the form of free nucleons n and p. For T > 1 
MeV and thus t < 1 sec, weak interactions are rapid (rates per nucleon r„^. p S> 
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H) and thus 



nv e -H> pe (4) 

pv e o ne + (5) 



drive neutron and protons to an equilibrium ratio 



— = e -( m n-m p )/l _ /g-j 

Thinking of the nucleon as a two-level system, this simply the Boltzmann ratio 
of the excited to ground state populations. 

At T = Tf ~ 1 MeV, the n — p interconversion (eq. 0|) stop as the weak inter- 
action "freezes out" (T n ^ p <C H), fixing n/p « e ~( m ™~ m p)/ T t ^ 1/6. Deuterium 
production p(n, 7)d occurs, but is stymied by the large number of photons per 
baryon rij/n^ = I/77 ~ 10 9 , which leads to effective deuteron photodissociation 
by the E 1 > B& = 2.22 MeV tail of the Planck distribution. During this time, 
free neutron decay reduces the neutron-to-proton ratio to n/p ~ 1/7. 

At T ~ 0.07 MeV, blackbody photons become ineffective to photodissociate 
deuterium. The deuteron abundance rapidly rises, and from this all of the light 
elements are built via strong (i.e., nuclear) interactions. A simplified reaction 
network appears in Fig. [TJ highlighting the reactions which dominate production 
of the light nuclides. In contrast to much of stellar nucleosynthesis, for BBN 
the number of key reactions is small and well-defined, and all of the important 
reactions have been measured in the laboratory at the relevant energies; no low- 
energy extrapolations are needed. 

Figure [2] shows the standard BBN light-element abundances as a function of the 
single free parameter 7/10 = 10 10 77 (eq. [3]). The vertical yellow band is the WMAP 
rj range (see §2.3|) . We see that the 4 He abundance is weakly sensitive to 77 (note 
that the zero is suppressed in the top-panel abscissa). In contrast, deuterium 
drops strongly with 77 and 3 He decreases substantially. The 7 Li abundance is 
plotted after 7 Be decay and thus sums both mass-7 species. 7 Li production 
dominates the mass-7 abundance in the low rj regime of the plot, while 7 Be 
production dominates at the high-77 regime, leading to the 'dip" behavior. 

The envelopes around the curves in Fig. [2] correspond to the la uncertainties 
in the abundance predictions. These uncertainties are propagated from the error 
budgets-statistical and systematic-of the twelve dominant reactions shown in 
Fig. H] dHl[l9l[2niEIl[22l[23l[2ll[25l[26l[l2]). The uncertainties in 4 He are tiny 
(< 1%), those in D and 3 He are small (~ 7%), while the 7 Li uncertainties are the 
largest (~ 12% in the high-?/ regime of interest). 

Several aspects of lithium production are noteworthy. Mass-7 is produced 

both as 7 Li and as Be; the 7 Be — > 7 Li electron capture occurs long after BBN 
ceases (J27H . The 3 He(a,7) 7 Be channel dominates Be production. Destruction 
occurs via 7 Be(n,p) Li followed by the rapid 7 Li(p, a) 4 He reaction. Finally, 6 Li 
production in standard BBN is very small: 6 Li/H ~ 10 -14 , or 6 Li/ 7 Li < 10~ 4 
(E81E9D. 



2.2 Light Element Observations 

Measuring the primordial abundance of any light element remains challenging. 
The BBN levels set at z ~ 10 10 are reliably accessible only in sites at z < 3 and 
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oftentimes z ~ 0. Other nucleosynthesis processes have intervened, as evidenced 
by the nonzero metallicity of all measured astrophysical systems. Thus one seeks 
to measure light elements in the most metal-poor systems, and then to obtain 
primordial abundances requires extrapolation to zero metallicity . Our discussion 
will follow closely recent treatments in refs. (130( 131 ^132]) . 

2.2.1 Deuterium, Helium-3, Helium-4 Deuterium can be measured di- 
rectly at high redshift. It is present in distant neutral hydrogen gas clouds which 
are seen in absorption along sightlines to distant quasars. At present there are 
seven systems with robust deuterium measurements ([33 . 34 . 35 . 36 .l3T ll38|) . These 
lie around redshift z ~ 3 and have metallicity ~ 10~ 2 of solar; thus deuterium 
should be primordial. For these systems, 

— = (2.82 ±0.21) x 10~ 5 (7) 

H 

where the error has been inflated by a the reduced xl = 2.95. 

Helium-4 can be measured in emission from nearby metal-poor galaxies (extra- 
galactic H II regions). The challenge is to reliably infer abundances at the needed 
<J 1% level. Several recent analyses differ due to systematics in the extraction of 
abundances from nebular lines (139,40,41,42,431). The mass fraction of (1401) 



Y p = 0.249 ± 0.009 (8) 

has the largest and most conservative measure of the error budget; the allowed 
range overlaps with analyses of other groups. 

Helium-3 is at present only accessible in our Galaxy's interstellar medium (I44p . 
This unfortunately means it cannot be measured at low metallicity, and so its 
primordial abundance cannot be determined reliably (|45p : we will not use 3 He to 
constrain BBN. 

2.2.2 Lithium- 7 Lithium is measured in the atmospheres of metal-poor 
stars in the stellar halo (Population II) of our Galaxy. Due to convective motions, 
surface material in stars can be dragged to the hot interior where lithium is 
burned readily; this effect is seen in low lithium abundances in cool halo stars. 
Fortunately, the hottest (most massive) halo stars have thin convection zones, and 
show no correlation between lithium and temperature. We consider only lithium 
abundances in these stars. 

Figure [3] shows lithium and iron abundances for a sample of halo stars (|46p . 
Li/H is nearly independent of Fe/H; this flat trend is known as the "Spite plateau" 
in honor of its discoverers (|47p . But heavy elements such as iron ("metals") 
increase with time as Galactic nucleosynthesis proceeds and matter cycles in and 
out of stars. Thus the Spite plateau indicates that most halo star lithium is 
uncorrelated with Galactic nulceosynthesis and hence, lithium is primordial. 

Moreover, the Spite plateau level measures the primordial abundance. Thanks 
to a sustained effort of several groups (0Sll3SlH2llSai5Tll521l53ll53ll53l5Sll, a large 
sample of halo stars have measured lithium abundances. The dominant error are 
systematic. A careful attempt to account for the full lithium error budget found 

dSZD 

| = (1.23i°!) x ID" 10 (9) 

with this 95%CL error budget dominated by systematics (see also £ )3.ip . 



777 - °- 05 ( 10 ) 
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Finally, it is encouraging to note that lithium has now been seen in stars in 
an accreted metal-poor dwarf galaxy. The Li/H abundances are consistent with 
Spite plateau, pointing to its universality ([58|) . 

2.2.3 Lithium-6 Due to the isotope shift in atomic lines, 6 Li and 7 Li are in 
principle distinguishable in spectra. In practice, the isotopic splitting is several 
times smaller than the thermal broadening of stellar lithium lines. Nevertheless, 
the isotopic abundance remains encoded in the detailed shape of the lithium 
absorption profile. 

High-spectral-resolution lithium measurements in halo stars attain the preci- 
sion needed to see isotope signatures. 6 Li detections have been claimed, in the 
range (J4~6|) 

6 Li 
7 Li 

Fig. [3] shows the inferred 6 Li/H abundance for some of these stars; its constancy 
with metallicity is strikingly reminiscent of the ordinary Spite plateau and simi- 
larly seems to suggest a primordial origin. 

Lithium-6 observations remain controversial. It has been argued that stellar 
convective motions can alter the delicate lineshapes and mimic 6 Li (I59p . Thus 
there are only a few halos stars for which there is widespread agreement that 
6 Li has even been detected. Thus, the conservative approach is to take the 6 Li 
observations as upper limits, though it is of interest to see what is required to 
explain the " 6 Li plateau," if it exists. Regardlessly, the isotopic searches confirm 
that most of primordial lithium is indeed Li. 

2.3 Microwave Background Anisotropics as a Cosmic Baryome- 
ter 

It is difficult to overstate the cosmological impact of the stunningly precise CMB 
measurements by WMAP and other experiments. The temperature and polar- 
ization anisotropies encode a wealth of cosmological information. Temperature 
fluctuations robustly record acoustic oscillations of the (re) combining baryon- 
photon plasma within dark matter potential; for a review, see (|60p . One of the 
most precise and robust results is the measurement of the cosmic baryon density 
and thus of rj. 

Figure H] shows the sensitivity of the temperature anisotropy to the baryon 
density, as a function of angular scale (multipole) on the sky. Broadly speaking, 
increasing baryon density amplifies the odd peaks and depresses the even peaks. 
Accurate measurements of these peaks by WMAP and other experiments pins 
down the baryon density. The most recent 7-year WMAP data release gives 

rj = (6.19 ±0.15) x 10 -10 (11) 

a 2.4% measurement! 

2.4 Assessing Standard BBN: the Lithium Problem(s) Revealed 

Prior to WMAP, BBN was the premier means of determining the cosmic baryon 
density. Standard BBN has one free parameter, rj, but three light elements have 
well-measured primordial abundances: D, 4 He, and 7 Li. Thus the problem is 
overdetermined: each element ideally selects a given value of rj, but allowing for 
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uncertainties actually selects a range of rj. If the different ranges are concordant, 
then BBN and cosmology are judged successful, and the cosmic baryon density is 
measured. This method typically specifies r\ to within about a factor ~ 2 (|6ip62p . 

The exquisite precision of the CMB-based cosmic baryon density suggests a 
new way of assessing BBN (f9l fl0jl . We exploit the CMB precision by using r] mwap 
as an input to BBN. This removes the only free parameter in the standard theory. 
Propagating errors, we compute likelihoods for all of the light elements. Fig. [5] 
shows these likelihoods (|30p based on WMAP data (63). Also shown are measured 
primordial abundances as discussed above. 

Figure [5] shows that deuterium observations are in spectacular agreement with 
predictions-the likelihoods literally fall on top of each other. This concordance 
links z ~ 3 abundance measurement with z ~ 10 10 theory and z ~ 1000 CMB 
data, and represents a triumph of the hot big bang cosmology. We also see that 
4 He predictions are in good agreement with observations. And, as noted in §2.21 
no reliable primordial 3 He measurements exist. 

Turning to 7 Li, the BBN+WMAP predictions and the measured primordial 
abundance completely disagree: the predictions are substantially higher than the 
observations. Depending on the treatment of systematic errors in the measured 
Li/H, the discrepancy is a factor Libbn+wmap/Li bs = 2.4 — 4.3, representing a 
4.2 — 5.3(7 discrepancy. This substantial mismatch constitutes the lithium problem 
(i.e., the 7 Li problem). 

Finally, as noted in £j2j standard BBN predicts an unobservable 6 Li/H abun- 
dance and 6 Li/ Li ratio far below the putative 6 Li plateau. To the extent that the 
6 Li plateau is real, this would constitute a second Li problem-the 6 Li problem. 

Hereafter we will focus largely on the well-established Li problem, but where 
appropriate we will mention the 6 Li problem. 

3 Solutions to the Lithium Problem(s) 

We have seen that the lithium problem was brought into sharp relief with the 
advent of the WMAP era, has become increasingly acute since. Possible solutions 
fall in three broad classes corresponding which part of the preceding analysis is 
called into question: 

1. Astrophysical solutions revise the measured primordial lithium abundance. 

2. Nuclear Physics solutions alter the reaction flow into and out of mass-7. 

3. Solutions beyond the Standard Model invoke new particle physics or non- 
standard cosmological physics. 

We consider each in turn. 

3.1 Astrophysical Solutions 

We first consider the possibility that BBN predictions are sound, i.e., standard 
cosmology and particle physics are correct, and the nuclear physics of mass-7 
production is properly calculated. If so, then the measured value of the primordial 
lithium abundance must be in error. In particular, the true value must be higher 
by a factor 3-4. 

As described in §2.2|, lithium abundances are measured via absorption lines in 
the photospheres of primitive, low-metallicity stars. For each star, the lithium 
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line strength is used to infer the Li/H abundance. Lithium abundances are nearly 
insensitive to metallicity (Fig. [3])— this Spite plateau implies that lithium is inde- 
pendent of Galactic nucleosynthesis and is primordial, and the plateau level is 
taken as the primordial abundance. If missteps exist in this chain of reasoning, 
the lithium problem could potentially be alleviated. 

For example, systematic errors could shift Li/H ratios in each star. We seek the 
total lithium content of the stellar photosphere, i.e., summed over all ionization 
states. However, the single (!) accessible 670.8 nm lithium line is sensitive only to 
neutral Li°. But in the stars of interest, lithium is mostly singly ionized. One must 
therefore introduce a large ionization correction Li + /Li° that is exponentially 
sensitive to the stellar temperature. Thus, a systematic shift upward in the 
temperature scale for halo stars would increase all stellar Li abundances and 
raise the Spite plateau towards the BBN+WMAP prediction. 

In practice, accurate determination of stellar temperatures remains non-trivial, 
because the emergent radiation is not a perfect Planck curve (else no lithium 
lines!), nor is local thermodynamic equilibrium completely attained in the stellar 
atmospheres. Quantitatively, the needed systematic shift in the temperature 
scale is about AT e g ~ 500 — 600 K upward, a ~ 10% increase over fiducial values 
(outside of previous claimed errorbars) (|64p . A re-evaluation of one method to 
determine stellar temperatures indeed corrects the scale, typically by ~ +200 
K ([65]) . However, later detailed studies of the stellar temperature scale are in 
good agreement with the fiducial temperature scale, leaving the lithium problem 
unresolved ()54p 66p. 

An entirely separate question remains as to whether a star's present lithium 
content reflects its initial abundance in the star. If the halo stars have destroyed 
some of their lithium, their present Li/H ratio sets a lower limit to the primordial 
lithium abundance. Indeed, given the low nuclear binding of 7 Li, it need only 
be exposed to relatively low stellar temperatures (T > 2.5 x 10 6 K) to suffer 
substantial destruction over the many-Gyr lifespan of a halo star. 

Lithium depletion is a major diagnostic of stellar structure and evolution 
([67l[68|[69l[70ll . For stars of solar composition, lithium destruction ([7"T|l72~l) ([72]) 
has long been studied in stellar evolution models. The major effect is convection, 
which circulates photospheric material deep into the interior (though still far 
from the stellar core) where nuclear burning can occur. Models for the evolution 
of low-metallicity stars, appropriate for the Spite plateau, now include numerous 
mixing effects which can change the photospheric lithium: convective motions, 
turbulence, rotational circulation, diffusion and gravitational settling, and inter- 
nal gravity waves ([701 175 1 174 ^ F75j) . These effects must occur at some level, and 
models have some success in fitting some observed trends in halo stars. There 
is general agreement that for stars with low metallicities, convective zones are 
substantially shallower than in solar-metallicity stars, and so depletion such be 
much smaller than the factor ~ 10 2 in the Sun (|76p . However, it remains difficult 
for models to quantitatively fit all of the data (|49p . 

Thus observational efforts to find clues for lithium depletion in the Spite plateau 
stellar data themselves remain of utmost importance. One study found the Spite 
plateau in field halo stars to be very thin, with no detectable star-to-star varia- 
tions around the Li-Fe trend, which showed a small positive slope in Li/H vs Fe/H 
([48p . A small lithium increase with metallicity is required due to contamination 
from Galactic cosmic-ray production of 7 Li and 6 Li f[57|) . An analysis of lithium 
and iron abundances in stars from the same globular cluster found trends con- 
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sistent with lithium depletion via diffusion and turbulent mixing; some models 
suggest these effects could remove the lithium problem entirely (|50p . However, 
systematic differences between globular cluster and field star lithium abundances 
raise concern about globular clusters at sites for constraining primordial lithium 

Recently, several groups have found that at very low metallicity, [Fe/H] <s —3, 
lithium abundances on average fall below the Spite plateau (i.e., below the levels 
seen at metallicities — 3 <> [Fe/H] <^ —1.5) (|52p 53 ( l5"T jl55"j) . These groups also find 
that the star-to-star scatter in Li/H becomes significant below [Fe/H] < —3. This 
appears to confirm the presence of significant lithium depletion in at least some 
halo stars. 

The recent evidence for lithium depletion at very low metallicities is a major 
development, yet its implications for primordial lithium remain unclear. No sig- 
nificant scatter is detected in plateau stars with —3 < [Fe/H] < —1.5. Also, 
in no stars is Li/H seen above the plateau, and in no metal-poor stars is Li/H 
near the WMAP+BBN value. Finally, while 6 Li measurements are difficult and 
controversial, there is general agreement that 6 Li is present in at least some 
plateau stars. This much more fragile isotope strongly constrains thermonuclear 
burning processes-if stellar material is exposed to temperatures hot enough to 
significantly reduce 7 Li, 6 Li should be completely destroyed. (|77|) . 

To summarize, determination of the primordial lithium abundance continues 
to be the focus of rapid progress. At present, however, the observational status 
of primordial lithium remains unsettled. A purely astrophysical solution to the 
lithium problem remains possible. On the other hand, the observed lithium 
trends-particularly the small lithium scatter in temperature and metallicity, and 
the presence of 6 Li~strongly constrain (but do not rule out) solutions of this 
kind. Consequently, it is entirely possible that the lithium problem cannot be 
resolved astrophysically, and thus we are driven to seek other explanations of the 
discrepancy; we now turn to these. 

3.2 Nuclear Physics Solutions 

We now consider the possibility that the measured primordial lithium abundance 
is correct, and the Standard Model of particle physics and the standard cosmology 
are also sound. In this case, the lithium problem must point to errors in the BBN 
light element predictions, in the form of incorrect implementation of standard 
cosmological and/or Standard Model physics. 

However, the standard BBN calculation is very robust and thus difficult to 
perturb. As summarized in Sj2j standard BBN rests on very well-determined 
physics applied in a very simple system. The cosmological framework of BBN is 
that of a very homogeneous universe (guaranteed by the smallness of the observed 
CMB temperature fluctuations (j3j)), with a cosmic expansion governed by exact 
expressions in General Relativity. The microphysics is that of the Standard 
Model, also very well-determined. The relativistic species, which comprise cosmic 
radiation that dominates the energy density, are very well thermalized and thus 
their properties are that of Bose-Einstein and Fermi-Dirac gasses, for which exact 
expressions are also available. 

The weak and strong (i.e., nuclear) interactions are also well-grounded in the- 
ory and calibrated empirically, but for BBN the needed physics is complicated 
(nuclear networks are large) and lies the farthest from first principles. Thus, 
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these are the only possible "weak links" in the standard BBN calculation, and it 
is here that solutions to the lithium problem have been sought. 

3.2.1 New and Revised Reactions One possibility is that weak and 
nuclear reactions in the BBN calculations are miscalculated due to reactions that 
are entirely missing, or that are included but whose rates are incorrect. But as 
described in ^2] and seen in Fig. [H only a relatively small number of reactions have 
been found to be important for producing the light elements, and all of these have 
been measured in the laboratory at the relevant energies. Their uncertainties have 
also been calculated and propagated through the BBN code, and are folded into 
the likelihoods appearing in Fig. [5] Moreover, BBN calculations use a much more 
extended reaction network than the simplified view of Fig. [H with all initial state 
pairings of A < 7 species present but most practice unimportant (fT8"| l?9 | [2"3" l [80 | 
81,82,83). Thus, to change the primordial lithium predictions requires surprises 
of some kind-either (i) the cross sections the known important reactions have 
uncertainties far beyond the quoted errors, or (ii) the cross sections for normally 
unimportant reactions have been vastly underestimated. 

For the important reactions seen in Fig. [TJ mass-7 production is dominated by 
the single reaction 3 He(a,7) Be. While the quoted error budget in the measured 
cross section is small, ~ 7% (|84p . absolute cross sections are difficult to measure. 
However, this reaction is also crucial in the production of solar neutrinos. To 
fix the cosmic lithium problem, the 3 He(a,7) 7 Be normalization would need to 
be low by a factor 3 — 4; if this were the case, the Be and 8 B solar neutrino 
fluxes would be lower by a similar factor. Thus we can view the sun as a reactor 
which probes the 3 He(a, 7) 7 Be rate, and the spectacular and precise agreement 
between solar neutrino predictions and observations becomes a measurement of 
the rate normalization which confirms the experimental results and removes this 
as a solution to the lithium problem (I85|) . 

Weak rates in BBN have received a great deal of attention over the years 
(JMlEBEHJESlEOlEIlEaESlEil). The basic nop interchange rates (eq. 0J) are 
most accurately normalized to the neutron lifetime. Corrections to the tree-level 
rates have < 1% effects on abundances, and thus are far too small to impact the 
lithium problem (|94p . 

Corrections to the standard thermonuclear rates have been considered as well. 
The effects of nonthermal daughter particles has been studied and found to be 
negligible (l95| i82 ) . Plasma effects, and electron Coulomb screening are also unim- 
portant (|96|) . 

Turning then to (normally) subdominant reactions, Angulo et al. (|97p noted 
that the (nonresonant) cross section for 7 ~Be(d, a)ap was poorly determined and 
could solve the lithium problem if it were a factor ~ 100 larger. They measured 
the cross section at BBN energies, but found values a factor ~ 10 smaller than 
had been used. 

Finally, the possibility of entirely new reactions has been recently studied by 
Boyd et al. (|82|) . These authors systematically considered a large set of reactions, 
some of which have been neglected in prior calculations. The focus of this study 
is almost exclusively on nonresonant reactions, with the result that even when 
allowing for extremely large systematic uncertainties in known cross sections, 
most new channels remain unimportant. The loophole to this analysis is the 
presence of new or poorly measured resonances. 
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3.2.2 Resonances Both standard and nonstandard reaction pathways to 
primordial mass-7 are firmly anchored to experimental data. The only remain- 
ing new nuclear physics can only intervene via resonances which have evaded 
experimental detection or whose effects have been underestimated. Cyburt and 
Pospelov (98) point out that the production of the known resonance 7 Be + d — > 
9 B (16.71 MeV) is poorly constrained experimentally; this reaction appears in 
Fig. [TJ Within current uncertainties, this resonance could promote the Be + d 
channel to become the dominant 7 Be destruction mode, and solving the lithium 
problem in an elegant manner. 

Generalizing the Cyburt and Pospelov (|98|) suggestion, ref. ([99]) searched the 
entire resonance solution space for BBN. These authors systematically considered 
all compound states created in mass-7 destruction, via all possible 2-body reac- 
tions of the form (n,p, d, £, 3 He, 4 He) + (Li, Be). Most possibilities were found 
to be unimportant. However, in addition to the 7 Be + d — >• 9 B (16.71 MeV) 
resonance, two other potentially important states were identified. The 7 Be + 1 — > 
10 B (18.80 MeV) resonance is known and within present uncertainties could be 
significant. On the other hand, there is little data on high-lying states of 10 C, but 
if a 10 C (15.0 MeV) exists and has J w = l - or 2~, this also could bring cosmic 
lithium into concordance if the reaction widths are large enough. This last possi- 
bility would be a homage to Fred Hoyle's celebrated prediction of the 12 C (7.65) 
state which solved the "carbon problem" of stellar nucleosynthesis (|100p . 

Fortunately all of these states are experimentally accessible. To identify or 
exclude them marks the endgame for a nuclear solution to the lithium problem. 

3.3 Solutions Beyond the Standard Model 

Finally, we turn to the most radical class of solutions the lithium problem. 
Namely, we assume that primordial lithium has been correctly measured, and 
that the nuclear physics of BBN has been calculated correctly and holds no sur- 
prises. In this case, we are forced to question the assumptions underlying the 
standard BBN calculation, i.e,. we must go beyond the Standard Model of par- 
ticle physics and/or the standard cosmology. For details beyond the overview 
below, see refs. (116 pi5p . 

3.3.1 Dark Matter Decay and Supersymmetry The existence of dark 
matter is now well-established, and its cosmic abundance has been inferred pre- 
cisely; for recent reviews, see refs. (|10Hll02"|ll03p . If dark matter takes the form 
of a relic particle created in the very early universe, it must be nonbaryonic. 
No Standard Model particles have the right properties, and thus dark matter 
demands physics beyond the Standard Model. Dark matter must of course be 
present during BBN, but ordinarily is assumed to be both nonrelativistic and 
weakly interacting, hence unimportant to cosmic dynamics and microphysical 
interactions. Similarly, dark energy is assumed to be negligible. 

While the identity of the dark matter is unknown, a simple, popular, and 
physically well-motivated possibility is that dark matter today consists of relic 
weakly interacting massive particles (WIMPs). In particular, if the universe be- 
gins with equal abundances of WIMPs and anti- WIMPs (if the two are distinct), 
then their abundance today is determined by the freezeout of their annihilations. 
Famously, to reproduce f2 m ~ 0.3 today, the annihilations must occur at ~ TeV 
scales (well before BBN). By happy coincidence, this is also the scale of the weak 
interaction, and of current accelerator experiments-this is the 'WIMP miracle" 
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PU1U1U21[T03]) . 

Moreover, it is likely that WIMPs today are the stable endpoints of a de- 
cay cascade. If so, then the WIMPs are the daughters born in the decays of 
the next-lightest particles in the cascade. The nature of these decays is model- 
dependent, but in general produce Standard Model particles which interact with 
the background plasma. If the decays occur during or after BBN, the interac- 
tions can change light element abundances (|104| ll05 |Hn51Hn71HMIT09H110| HTT | 
I112pil3pil4llll5pil6llll7|) . and thus potentially could solve the lithium problem(s) 

([IiainillMlmllmllmlllMlElI) ■ 

We thus consider the effects of decaying massive particles during or after BBN. 
To get a feel for the basic physics, consider a particle X (and X if they are 
distinct) with mass mx S> Tn p , which decays with lifetime tx- The decays can 
have electromagnetic and/or hadronic channels, and given the massive nature of 
X, these decay products will be very relativistic and thus nonthermal. 

As these electromagnetic and/or hadronic cascades thermalize in the cosmic 
environment, they interact with the light elements, largely via fragmentation 
(photodissociation or spallation). For example, high energy nucleons N fragment 

4 He: 

2n + 2p + N 

This reduces the 4 He abundance, but more importantly creates new deuterium 
and 3 He. Furthermore, the secondary particles in eq. f)12|) are themselves non- 
thermal, and can initiate further interactions with the background thermal light 
nuclides. Of particular importance is the conversion of 7 Be into 7 Li with sec- 
ondary neutrons 

n + 7 Be^p + 7 Li , (13) 

which substantially enhances mass-7 destruction because of the lower Coulomb 
barrier for Li. Also of interest is the nonthermal production of 6 Li via secondary 
nonthermal deuterons 

d + 4 He^ 6 Li + 7 . (14) 

Clearly, both of these processes are of great interest for both lithium problems. 
Figure [6] shows light-element abundances in the presence of a hadronically 
decaying particle X. Abundance contours are plotted as a function of decay 
lifetime tx, and the pre-decay X abundance 

m x n x n x n -, 

Cx = = rn x — V (15) 

n 7 rtb 

In the absence of X decays, cosmic expansion dilutes baryons and X particles in 
the same way, and thus Cx remains constant. Figure [6] is for a hadronic branching 
fraction B^ = 1, i.e., all decays produce hadrons. 

In each panel of Fig. [6l the colored areas indicate parameter regions where 
the predicted light element abundances disagree with observations, while the 
remaining white regions are allowed. At fixed tx, the light element perturbation 
is proportional to the X abundance. The limit Qx — > lies at the bottom of 
the plot, and corresponds to the unperturbed (standard BBN) case. At fixed 
abundance Cx, the light element effects strongly depend on particle lifetime, 
which determines when the light elements are perturbed. For example, at tx ^ 10 
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sec, decays occur before the light elements are formed, and so the light elements 
are unaffected. Hadronic decays dominate the perturbations in the tx ^ 10 6 sec 
of most interest here; electromagnetic decays dominate at longer times. 

In the 7 Li/H panel, the colored regions include the unperturbed low-^x regime, 
reflecting the lithium problem for standard BBN. However, for a relatively narrow 
region with lifetimes Tx ~ 10 2 — 10 3 sec, that the 7 Li abundance is reduced 
and brought into accord with observation. This arises due to the production 
of secondary neutrons (eq. fT2l) . which facilitate mass- 7 destruction via 7 Be-to- 
7 Li conversion (eq. PT3]k at longer lifetimes the secondary neutrons decay. For 
10 2 — 10 3 sec lifetimes to be viable, however, the other light elements must remain 
in concordance. As seen in the final, summary panel, in the regime in which 7 Li 
is reduced, all other constraints are satisfied except D/H, which is unacceptably 
high due to secondary deuteron production. This tension between deuterium 
and 7 Li is a fundamental feature of decay scenarios, and as we see may allow for 
solutions but requires fine tuning. 

Supersymmetry provides well-motivated candidates for decaying dark matter 
(|125pi26]ll27p . Supersymmetry doubles the particle content of nature by re- 
quiring opposite-statistics (fermion o boson) partners for every known particle. 
These partners are produced copiously in the very early universe. The light- 
est supersymmetric partner (LSP) is the stable end product of the decays of 
higher-mass supersymmetric particles, and naturally becomes a dark matter can- 
didate (|128|I129| . These decays are a fundamental aspect of supersymmetric 
dark matter, and thus supersymmetric scenarios demand that particle decays 
occur. Moreover, even in the simplest scenarios-i.e., the constrained minimal 
supersymmetric Standard Model-the lifetime of the next-to-lightest supersym- 
metric partner can be long: > 100 sec. Thus, minimal supersymmetry holds the 
tantalizing possibility of solving the lithium problem. 

Figure [7] illustrates the interplay between supersymmetry and the lithium prob- 
lem (|32h . The context is minimal supersymmetry in which the spin-3/2 gravitino 
is the next-to-lightest partner, which decays into the LSP neutralino that com- 
prises the dark matter today. For a benchmark scenario in this model, the su- 
perpartner mass spectrum, lifetimes, and decay products are calculated. These 
are used to compute light element abundances, with nuclear uncertainties prop- 
agated; the results are compared to observations. The resulting x 2 1S plotted as 
a function of gravitino mass m^/? and pre-decay abundance £3/2 ■ At low C3 1^ , we 
recover the unperturbed, standard case, where x 2 ~ 32 for one effective degree 
of freedom (3 light elements — 2 parameters), a poor fit. At high £3/21 the light 
element perturbations are worsened over the standard case, and these supersym- 
metric parameters are excluded by BBN. Most interestingly, along the diagonal 
loop, the fit improves, and in the interior "islands," x 2 drops below 6, correspond- 
ing to ~ 2.4o". This region shows substantial improvement over standard BBN, 
and physically arises as the regime of optimal tradeoff between 7 Li destruction 
and deuterium production. This region thus tantalizingly stands as a possible 
supersymmetric solution to the lithium problem, albeit statistically marginal and 
fine-tuned. 

Finally, in very recent years an entirely new aspect of decaying dark matter 
alteration of BBN has emerged. If the decaying particles are electrically charged, 
then the negatively charged dark matter can form bound states with charged 
nuclei, such as (pX~), ( 4 HeX~), and ( 7 BeX~). As pointed out by Pospelov 
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(|117p . these bounds states unleash a rich array of new effects, in addition to 
the perturbations accompanying X~ decays. For the heavy X~ of interest, the 
binding energy of (A z X~) is B = Z 2 Aa 2 m p /2 sa 30Z 2 A keV, comparable to nu- 
clear binding. The Bohr radius is a = (2aZAm p )~ 1 ~ 1A~ 1 Z~ 1 fm, comparable 
to nuclear sizes. Bound-state Coulomb barriers are thus reduced and ( 7 BeA"~) 
becomes easier to destroy. 

Beyond these basic effects, a new "bound state chemistry" can occur, whose 
nature is still under active study. One important effect (1117P is catalysis, par- 
ticularly d + ( 4 HeX) — > 6 Li + X which enhances 6 Li production far above the 
ordinarily small radiative capture d(a, 7) 6 Li. This can enhance 6 Li production by 
orders of magnitude, possibly addressing the 6 Li problem if real, but oftentimes 
overproducing 6 Li. Rates for catalyzed reactions have recently been calculated 
to high precision (|130p . Intriguingly, ( 8 BeX~) states have a binding energy very 
close to the energy for 8 Be — > aa breakup; if the binding energy is larger, then 
( 8 BeX~) is stable and can allow for cosmological production of 9 Be ([12211131"]) . 

The general properties of X~ recombination and bound states, and their impact 
on BBN, have been studied in refs. (fTTTlll32lll23ll33lll22l[15]) . Looking at bound 
state effect only (i.e., ignoring decay effects), catalyzed 7 Be destruction is effective 
and catalyzed 6 Li is not overproduced, for sufficiently large X~ abundance in 
the regime tx ~ 2000 sec f)132[) . Indeed, there exist regions of parameter space 
wherein both the Li and Li problems are solved. Full calculation of bound state 
and decay effects together is required to verify if solutions remain in specific 
detailed models; early calculations confirm solution space exists around tx ~ 10 3 
(H34J1123J. 

Catalyzed production of 9 Be also occurs and is constrained by the non-observation 
of a beryllium "Spite plateau;" the resulting limits on tx are comparable to those 
imposed by limits on primordial 6 Li ()122p . 9 Be constraints have the added advan- 
tage that they rely on elemental abundances which are much simpler to obtain 
reliably than isotopic abundances. However, very recently 9 Be was investigated 
using updated catalysis rates from ref. (|130p . which greatly reduce the 9 Be pro- 
duction rate. The resulting 9 Be abundance then is quite small, below observable 
levels (fT3T]) . 

Bound-state effects are important in the context of minimal supersymmetry. 
Substantial parameter space exists in which the gravitino is the LSP and thus the 
dark matter, while the next-to-lightest partner is the charged stau f , the scalar 
partner to the tau lepton. These models are probed by bound-state BBN, which 
imposes constraints that are complementary to accelerator limits (]121[I132"]I135[ 
[T33l[T34l[T22l[T24l[T36D . Indeed, the solutions to the 7 Li and possibly also the 6 Li 
problem which exist at Tx ~ 1000 sec could be interpreted as support for these 
models. 

To summarize, decaying dark matter scenarios introduce a rich array of novel 
processes that can alter light elements during and after BBN. Moreover, such 
scenarios find well-motivated origin in super symmetric cosmologies. Indeed, 
decaying-particle BBN offers important constraints on supersymmetry. Further- 
more, the Li and possibly also 6 Li problems can be solved in decaying particle 
scenarios which are realized in plausible minimal supersymmetric scenarios. This 
area is ripe for further theoretical, observational, and experimental development. 
3.3.2 Changing Fundamental Constants Observations of multiple atomic 
transitions in metals residing in high-redshift quasar absorption systems test fun- 
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damental physics in environments at great spacetime separations from our own; 
for a review see ref. (|137p . Surprisingly, some data hint at variations in the fine- 
structure constant at z ~ 3, showing Soem/oiem — —0.5 x 10 -5 at the ~ 5<r level, 
while others find results consistent with no variation. Thus the observational sit- 
uation is unsettled, but intriguing. 

Time variations in low-energy physics can be accommodated and are even 
expected in the context of some unified theories which in general predict stronger 
variations at earlier times. Moreover, an underlying unified theory implies that 
all Standard Model couplings and particle masses should vary, with definite but 
model-dependent interrelationships. 

There is thus theoretical impetus, and some observational motivation, to con- 
template changes in fundamental constants during BBN. The change in light 
elements depends on which parameters (couplings, masses) change, and on the 
size of the perturbations (1138, 139 lll40yi4iyi42p . In general, there is large model- 
dependence in quantifying the these variations, the links among them, and their 
manifestation in nuclear properties (masses, binding energies, cross sections). 

An alternative approach is to turn the problem around and to consider the BBN 
implications of variation in nuclear physics parameters (|143pi44| . Coc et al. (|140p 
systematically study the light-element response to variations in oem, the electron 
mass m e , the neutron lifetime r n , the neutron- proton mass difference m n — m p , 
and the deuteron binding energy Br). Of these, the most sensitive parameter is 
the deuteron binding energy. A change of —0.075 < SBd/Bo < —0.04 lowers the 
Li abundance into concordance, without perturbing 4 He or D/H beyond their 
observed error range. Thus, unified models which predict changes of this order 
can solve the lithium problem. 

3.3.3 Nonstandard Cosmologies The lithium problem could indicate 
nonstandard cosmology rather than particle physics. One recent such proposal 
is that cosmic acceleration could result from large-scale inhomogeneities in the 
cosmic density. Isotropy constraints can be satisfied if we occupy a privileged 
view from nearly the center of a spherically symmetric cosmic underdensity, which 
only returns to the cosmic mean at horizon-scale distances ()145|) . Such a scenario 
explains cosmic acceleration within General Relativity, and without invoking dark 
energy, but must abandon the cosmological principle, instead requiring that ours 
is a privileged view of the universe. 

Esthetics aside, observations probe such an inhomogeneous cosmos. BBN oc- 
curs differently in such a universe, if the baryon-to-photon ratio varies along the 
inhomogeneity. In particular, ref. f|146|) emphasizes that the observations of 7 Li 
are made locally, at low z, while D/H and the CMB are both measured in the 
distant universe at high z. If the local baryon-to-photon ratio 770 is low by a factor 
~ 2, then indeed one would expect local 7 Li to fall below the WMAP prediction, 
while D/H would agree. 

This clever scenario however must face an array of observational tests. Re- 
laxed galaxy clusters probe the cosmic baryon-to-matter fraction, and show vari- 
ations < 8% out to z ~ 1, far less than the ~ 50% variation needed for lithium 
abundances (I147p . Moreover, local measurements of D/H even more directly 
constrain the local 77 measurement. Because some stellar destruction may well 
have occurred, the local D/H sets a lower limit to the primordial abundance. In 
the lower halo of our own Galaxy, D/H = (2.31 ± 0.24) x 10 -5 . This value is 
in good agreement with the high-z D/H measurements and the WMAP+BBN 
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predictions; it is therefore inconsistent with the low r\ needed for a low local 
7 Li/H. 

Other nonstandard cosmology scenarios have been proposed to solve the lithium 
problem. One suggests that a large fraction (~ 1/3 — 1/2) of baryons were pro- 
cessed through the first generation of stars (Population III) which ejected lithium- 
free matter (|148p . Such models face grate difficulties due to the substantial D/H 
depletion and 3 He production which must also occur (|149pi50|) . 

4 Discussion and Outlook 

BBN has entered the age of precision cosmology. This transition has brought 
triumph in the spectacular agreement between high-redshift deuterium and the 
BBN+WMAP prediction, and in the WMAP confirmation of the longstanding 
BBN prediction of nonbaryonic dark matter. But new precision has raised new 
questions: the measured primordial 7 Li abundance falls persistently and signifi- 
cantly below BBN+WMAP predictions. Moreover, there are controversial hints 
of a primordial 6 Li abundance orders of magnitude above the standard prediction. 
As we have seen, disparate explanations for the lithium problem(s) remain 
viable. Fortunately, most alternatives are testable in the near future. 

1. Astronomical observations. Recent indications of lithium depletion in ex- 
tremely metal-poor halo stars are tantalizing. In the coming Great Survey 
era, we may expect many more such stars to be identified, and the lithium 
trends explored in large statistical samples. These will require careful com- 
parison with theory. Observations of 6 Li remain challenging, and as yet it 
remains unclear what trends exist with metallicity. 

Great insight would result from alternative measures of primordial lithium, 
e.g., in the interstellar medium of metal-poor galaxies nearby or at high 
redshift. 

2. Nuclear experiments. The enormous effort of the nuclear community has 
empirically pinned down nearly all nuclear inputs to BBN. Remaining are 
a few known or proposed resonances which would amplify 7 Be destruction. 
These are within reach of present facilities, so that the nuclear physics of 
standard BBN can and will be fully tested. 

3. Collider and dark matter experiments. The LHC is operational and much 
of minimal supersymmetry lies within its reach. The discovery of super- 
symmetry would revolutionize particle physics and cosmology, and would 
transform decaying particle BBN scenarios into canonical early universe 
cosmology. Alternatively, if the LHC fails to find supersymmetry and/or 
finds surprises of some other kind, this will represent a paradigm shift for 
all of particle physics and particle cosmology, and BBN will lie at the heart 
of this transformation. 
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Figure 1: Simplified BBN nuclear network: 12 normally important reactions 
shown in blue, and proposed/tested new reactions in red. 
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Figure 2: BBN theory predictions for light nuclide abundances vs baryon- 
to-photon ratio rj. Curve widths: lcr theoretical uncertainties. Vertical band: 
WMAP determination of rj. 
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Figure 3: Lithium abundances in selected metal-poor Galactic halo stars, from 
(j46[) with permission. Fore each star, elemental Li = 6 Li + 7 Li is plotted at the 
star's metallicity [Fe/H] = log lo [(Fe/iJ) obs /(Fe/H) ]. The flatness of Li vs Fe 
is the "Spite plateau" and indicates that the bulk of the lithium is unrelated 
to Galactic nucleosynthesis processes and thus is primordial. The horizontal 
band gives the BBN+WMAP prediction; the gap between this and the plateau 
illustrates the Li problem. Points below the plateau show 6 Li abundances; the 
apparent plateau constitutes the 6 Li problem. 
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Figure 4: CMB sensitivity to cosmic baryon content. Predictions for tempera- 
ture anisotropy (rms of temperature fluctuation Ay = (Tp) 2 — (T) 2 ) plotted as 
a function of angular scale (multipole £, roughly corresponding to angular size 
180°/£). Baryon density is seen to be encoded in the values and particularly the 
ratios of the peak heights. Figure from (|60p . with permission. 
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Figure 5: Comparison of BBN+WMAP predictions and observations, from fj30|) . 
Plotted are likelihood distributions for light nuclide abundances. Blue curves: 
theory likelihoods predicted for standard BBN using the cosmic baryon density 
determined by WMAP (|63p . Yellow curves: observational likelihoods based on 
primordial abundances as in §2.21 Dotted curves: observational likelihoods for 
different analyses of abundance data; the difference between these and the yellow 
curves gives a sense of the systematic errors. Note the spectacular agreement of 
D/H, and in contrast the strong mismatch between 7 Li theory and data, which 
constitutes the lithium problem. 
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Figure 6: Effect of nonthermal particle injection on light element abundances, 
from pip . Each panel shows abundance contours in the presence of the hadronic 
decay of a (neutral) particle X, plotted as a function of X abundance Qx ( e( i- fl5|) 
and mean life tx- Results are shown for a hadronic branching ratio Bh = 1. As 
summarized in the last panel, the parameter regions where the Li problem is 
solved also lead to deuterium production, placing the two in tension. 
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Figure 7: Light element constraints on gravitino decays in the context of the 
constrained minimal supersymmetric Standard Model, from (|32p . Shown are x 2 
contours based on light element abundance constraints and 1 effective degrees of 
freedom, over the space of gravitino mass m 3 / 2 an d abundance C3/2 (eg. [To]) . The 
cross indicates the maximum x 2 - diagonal "archipelago" shows the region where 
the lithium problem is greatly reduced by trading off 7 Li destruction for some 
degree of deuterium production. 



